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Preface 


This  Interagency  Report  (IR)  is  one  of  a series  of  IRs  that  will  form  the  basis  of  a method  for 
analysis  of  airborne  asbestos  by  transmission  electron  microscopy.  The  form  and  style  of  the 
American  Society  for  Testing  and  Materials  (ASTM)  was  adopted  as  a standard  format  for  this 
series  of  reports. 

Shirley  Turner  or  Eric  Steel  can  be  contacted  about  obtaining  the  programs  on  disk  or  from  a 
planned  world  wide  web  site. 


11 


1.  Scope 

1 . 1 This  test  method  describes  a procedure  for  determining  the  component  of  uncertainty  in 
the  concentration  of  asbestos  that  is  due  to  variation  in  sampling  of  a population  of  values 
(corresponding  to  the  Type  A evaluation  of  uncertainty  described  in  TN  1297^).  The  component 
of  uncertainty  arising  from  systematic  effects  (Type  B evaluation  of  uncertainty)  is  not  addressed 
in  this  method. 

1 .2  The  test  method  describes  use  of  a bootstrapping  procedure  for  determination  of  a 95% 
confidence  interval  for  the  concentration  of  asbestos  deposited  onto  filters. 

1.3  The  test  method  is  applied  where  the  concentration  of  asbestos  is  determined  by 
transmission  electron  microscopy. 

1.4  This  method  for  determination  of  a confidence  interval  is  nonparametric,  i.e.,  no 
probability  distribution  is  assumed  for  the  concentration  of  asbestos  in  air  or  on  filters. 

1.5  The  method  can  be  used  to  derive  95%  confidence  intervals  for  analyses  of  both  indoor 
and  outdoor  sets  of  filters  as  described  in  the  AHERA  method. 

1.6  The  method  can  be  used  to  derive  confidence  intervals  for  the  number  of  structures  per 
grid  square,  the  number  of  structures  per  mm^,  the  number  of  structures  per  cm^  of  air  or  for  any 
other  relevant  unit. 

1.7  The  method  cannot  be  applied  to  data  sets  in  which  all  values  are  zero. 

1.8  The  method  requires  analysis  of  a minimum  of  four  grid  squares  on  each  of  five  filters  and 
analysis  of  an  approximately  equal  number  of  grid  squares  per  filter. 

1.9  This  method  requires  use  of  an  Apple,  PC  (EBM  or  IBM  clone)  or  a computer  that  can 
compile  and  execute  the  code  in  Appendix  X4  and  X5. 

1.10  This  method,  though  written  specifically  for  asbestos,  can  be  applied  to  other  data  sets 
containing  at  least  20  values. 

2.  Terminology 

2.1  Definitions: 

2.1.1  bootstrapping— ?i  statistical  procedure  whereby  a set  of  data  is  uniformly  and  randomly 
sampled  with  replacement.  Statistical  estimates  are  then  determined  from  the  randomly  resampled 
data  sets. 

2. 1 .2  grid  square,  grid  opening— dj\  area  on  a grid  used  for  analysis  of  asbestos  by 
transmission  electron  microscopy. 

2.1.3  95%  confidence  interval— di  statistical  interval  derived  from  a random  sample  of  a 
population  that  may  enclose  a value  of  interest,  e.g.,  the  mean  of  a population.  If  95%  confidence 
intervals  are  repeatedly  calculated  from  many  independent  random  samples  of  a population,  the 
true  value  of  interest  will  be  contained  within  the  confidence  intervals  95%  of  the  time. 

2.2  Description  of  Terms  Specific  to  This  Standard: 

2.2.1  AHERA  method— mtihod.  for  analysis  of  asbestos  by  transmission  electron  microscopy 
developed  by  the  Environmental  Protection  Agency^  with  subsequent  modifications  by  the 
National  Institute  of  Standards  and  Technology.  This  method  was  developed  in  response  to  the 
Asbestos  Hazard  Emergency  Response  Act  (AHERA)  of  1986. 


^Taylor,  B.N.  and  C.E.  Kuyatt,  1994,  NIST  Technical  Note  1297,  1994. 
^Code  of  Federal  Regulations,  1987,  52,  No.  210,  41826-41905. 
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2.2.2  bootstrap  program— Xh&  computer  program  (for  either  the  Apple  or  PC)  associated  with 
this  test  method. 

2.2.3  input  data  set—z.  data  set  containing  the  information  to  be  input  into  the  bootstrap 
program.  The  data  set  consists  of  one  column  containing  an  identification  number  or  text  for  each 
grid  opening  analyzed  and  a second  column  containing  the  number  (or  concentration)  of 
structures  found  on  each  grid  opening. 

2.2.4  bootstrap  data  set—z  basic  set  of  values  chosen  by  uniform  random  resampling  in  the 
bootstrap  program.  In  the  bootstrap  program,  the  size  of  the  bootstrap  data  set  is  equal  to  the 
number  of  grid  squares  analyzed. 

2.2.4. 1 Discussion— bootstrap  program  determines  the  average  number  of  structures  per 
grid  square  for  each  bootstrap  data  set  and  stores  this  value  in  an  internal  data  file. 

2.2.5  number  of  iterations— Xh&  number  of  bootstrap  data  sets  generated  by  the  program. 

2.2.5. 1 Discussion— Th&  number  of  iterations  is  set  at  6000  for  this  program. 

2.2.6  grand  mean—Xht  average  value  of  the  mean  values  generated  for  the  bootstrap  data  sets 
from  a run. 

2.2.7  lower  95%  confidence  limit— Xht  lower  value  of  the  95%  confidence  interval  for  the 
average  number  of  structures  generated  for  the  bootstrap  data  sets  from  a run. 

2.2.7. 1 Discussion— In  this  program,  the  mean  values  generated  for  6000  bootstrap  data  sets 
are  sorted.  The  lower  confidence  limit  is  the  mean  value  that  is  150  values  from  the  lowest  value 
in  the  list. 

2.2.8  upper  95%  confidence  limit— X\iQ  upper  value  of  the  95%  confidence  interval  for  the 
average  number  of  structures  generated  for  the  bootstrap  data  sets  from  a run. 

2.2.8. 1 Discussion— In  this  program,  the  mean  values  generated  for  6000  bootstrap  data  sets 
are  sorted.  The  upper  confidence  limit  is  the  mean  value  that  is  150  values  from  the  highest  value 
in  the  list. 

2.2.9  run—zn  execution  of  the  program. 

2.2.9. 1 Discussion— Tht  output  from  each  run  consists  of  the  average  number  of  structures  in 
the  input  data  set,  and  the  grand  mean,  and  lower  and  upper  values  of  the  average  number  of 
structures  derived  from  the  bootstrap  data  sets. 

3.  Significance  and  Use 

3.1  The  uncertainty  derived  by  this  method  is  useful  for  a general  comparison  of  analyses 
obtained  by  different  operators  or  laboratories  on  the  same  sampling  area. 

3.2  The  uncertainty  derived  by  this  method  provides  an  indication  of  how  an  analytical  value 
compares  with  values  set  by  government  regulations. 

3.3  The  test  method  should  be  used  for  analytical  values  derived  from  the  AHERA  method. 
Where  all  values  are  zero,  the  confidence  interval  should  be  derived  assuming  a Poisson 
distribution. 

3.4  The  uncertainty  generated  by  this  method  should  be  included  in  reports  of  analysis 
resulting  from  use  of  the  AHERA  method. 

4.  Procedure 

4. 1 Obtain  an  analysis  of  air  samples  using  transmission  electron  microscopy  as  described  in 
the  AHERA  method. 
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4.2  Create  a file  containing  an  input  data  set  (see  example  input  data  sets  for  Apple  computers 
in  Figure  XI . 1 of  Appendix  XI  and  for  the  PC  in  Figure  X2. 1 of  Appendix  X2).  Two  separate 
files  should  be  created  and  run  separately  if  there  are  both  an  inside  and  outside  set  of  filters. 

NOTE  1— The  input  file  for  the  programs  is  an  ASCII  file. 

4.2.1  The  first  three  lines  of  the  input  data  set  contain  header  information  of  the  laboratory's 
choosing. 

4.2.2  The  fourth  line  contains  the  units  of  the  input  data.  Examples  of  units  include: 
structures,  structures/mm^,  structures/cm^,  etc. 

4.2.3  The  input  for  the  data  starts  on  the  fifth  line.  There  are  two  columns  for  the  data  set. 

The  first  column  contains  identification  numbers  or  text  for  the  grid  openings  analyzed  and  the 
second  column  contains  the  number  or  concentration  of  asbestos  structures  identified  in  each  grid 
opening. 

NOTE  2~The  entries  for  the  columns  cannot  exceed  8 characters.  The  columns  can  be  space  or 
tab  delimited. 

4.3  Run  the  bootstrap  program.  For  information  purposes,  a general  outline  of  the  program  is 
given  in  Appendix  X3  and  the  source  codes  are  given  in  Appendix  X4  and  Appendix  X5.  The 
Apple  version  is  run  by  double  clicking  on  the  program  icon  or  name  and  the  PC  version  is  run  by 
typing  the  program  name  at  the  DOS  prompt. 

4.3.1  Input  the  file  name  for  the  input  data  set  (see  Figure  XI .2  and  Figure  X2.2). 

4.3.2  Wait  while  the  program  is  performing  the  bootstrapping  procedure  until  the  statistics  for 
the  input  and  bootstrap  data  sets  are  displayed.  The  time  required  for  the  calculation  varies 
depending  on  the  type  and  speed  of  the  computer  used. 

4.3.3  Respond  to  questions  concerning  saving  the  output  data  to  a file  (see  Figures  XI .3  and 
X2.3  for  example  output  files). 

4.4  Use  a statement  identical  or  similar  to  that  given  in  the  output  data  file  for  reporting  the 
concentration  determined  using  this  method. 

NOTE  3 “The  number  of  significant  figures  given  by  the  bootstrapping  program  is  larger  than 
appropriate  for  most  analyses.  Laboratory  personnel  should  determine  the  number  of  significant 
figures  appropriate  to  report  for  the  analysis. 

5.  Precision  and  Bias 

5. 1 Precision--T>2iidi  has  been  collected  only  on  the  repeatability  of  determination  of  the  95% 
confidence  interval.  Data  on  reproducibilty  will  be  collected  later. 

5.1.1  Repeatability— Tht  repeatability  of  this  test  method  is  dependent  upon  the  nature  of  the 
data  set.  The  variation  in  the  confidence  limits  obtained  by  repeating  the  bootstrap  procedure  50 
times  is  given  for  two  examples  below. 

5. 1 . 1 . 1 Confidence  intervals  were  determined  for  the  data  set  given  in  Appendix  XI . For  the 
lower  95%  confidence  limit,  a value  of  6.7  s/mm^  was  obtained  for  31  of  50  bootstrap  runs  and  a 
value  of  10  s/mm^  was  obtained  for  the  other  19  runs.  For  the  upper  95%  confidence  limit,  a 
value  of  46.7  s/mm^  was  obtained  3 times  and  a value  of  50.0  s/mm^  was  obtained  47  times. 
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5. 1 . 1 .2  Confidence  intervals  were  determined  for  the  data  set  given  in  Appendix  X2.  For  the 
lower  95%  confidence  limit,  values  of  2.34,  2.36  and  2.38  st./grid  square  were  generated  15,  34 
and  1 time,  respectively.  For  the  upper  95%  confidence  limit,  values  of  3.40,  3.42  and  3.44 

St. /grid  square  were  generated  21,  28  and  1 time,  respectively. 

5.2  Bias—lhQ  bias  in  the  determination  of  confidence  intervals  from  bootstrapping  methods  is 
discussed  in  Efron  and  Tibshirani  (1993). 

6.  Keywords 

6. 1 asbestos;  bootstrapping;  transmission  electron  microscopy;  confidence  interval 


Additional  References 

Efron,  B.E.  and  R.J.  Tibshirani  (1993)  An  Introduction  to  the  Bootstrap.  Chapman  and  Hall,  NY, 
NY. 

Hahn,  G.J.  and  W.Q.  Meeker  (1991)  Statistical  Intervals:  A Guide  for  Practitioners.  Wiley- 
Interscience,  NY. 
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APPENDIX 


XI.  BOOTSTRAPPING  PROGRAM  RUN  ON  AN  APPLE  COMPUTER 
Fig.  XI.  1 Input  file  for  the  bootstrapping  program. 

This  is  a test  file  for  the  bootstrap  program. 


structures/mm2 
1 0 
1 0 
1 0 
1 100 
1 0 
1 0 
2 0 
2 0 
2 0 
2 0 
2 0 
2 0 
3 0 
3 0 
3 0 
3 0 
3 0 

3 0 

4 0 
4 0 

4 200 
4 0 
4 100 

4 100 

5 0 
5 0 
5 0 

5 100 
5 0 
5 200 
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Fig.  XI.  2 On  screen  input  and  output  obtained  during  a run  of  the  Apple  version  of  the 
bootstrapping  program. 


NIST  data  analysis  program. 

Hit  Return  and  select  a file  to  open. 

Calculating  bootstrap  - 


For  input  data  set: 


Mean  value 

= 

26.6667 

bootstrap  data  sets: 

Grand  mean  value 

= 

26.5894 

Lower  95%  confidence 

limit  = 

6.6667 

l^per  95%  confidence 

limit  = 

50.0000 

The  following  would  be  an  appropriate  way  to  state  the  results  of 
this  calculation: 

The  mean  value  for  this  data  set  is  26.6667  structures/irim2 . The  interval 
from  6.6667  to  50.0000  st2ructures/mm2  is  a 95%  confidence  interval  for  the 
mean  value;  i.e.,  we  may  assert  with  95%  confidence  that  the  mean  value 
lies  between  6.6667  and  50.0000  struct\ires/mm2 . 

Hit  Return  and  select  a file  for  output  or  cancel. 


6 


Fig.  XI  .3  Output  file  for  the  Apple  version  of  the  bootstrapping  program. 


example.out  - 6/13/95 

Input  data  set  = example 
Number  of  grid  openings  = 30 

This  is  a test  file  for  the  bootstrap  program. 


structures/mm2 

For  input  data  set: 

Mean  value  = 26.6667 

For  bootstrap  data  sets: 

Grand  mean  value  = 26.5894 

Lower  95%  confidence  limit  = 6.6667 

Upper  95%  confidence  limit  = 50.0000 

The  following  would  be  an  appropriate  way  to  state  the  results  of 
this  calculation: 

The  mean  value  for  this  data  set  is  26.6667  structures/mm2.  The  interval 
from  6.6667  to  50.0000  structures/mm2  is  a 95%  confidence 
interval  for  the  mean  value;  i.e.,  we  may  assert  with  95%  confidence 
that  the  mean  value  lies  between  6.6667  and  50.0000 
structures/mm2. 
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X2.  BOOTSTRAPPING  PROGRAM  RUN  ON  A PC 


Fig.  X2. 1 Input  file  for  the  bootstrapping  program. 


Example  data  set 


st./grid  square 
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Fig.  X2.2  On  screen  input  and  output  obtained  during  a run  (PC  version) 


NIST  data  analysis  program 
Input  the  file  name:  example3.dat 
Calculating  bootstrap  - 
For  input  data  set: 

Mean  value  = 2.8800 

For  bootstrap  data  sets: 

Grand  mean  value  = 2.8765 

Lower  95%  confidence  limit  = 2.3600 

Upper  95%  confidence  limit  = 3.4000 


The  following  would  be  an  appropriate  way  to  state  the  results  of 
this  calculation: 

The  mean  value  for  this  data  set  is  2.8800  st./grid  square.  The  interval 
from  2.3600  to  3.4000  st./grid  square  is  a 95%  confidence  interval  for  the 
mean  value;  i.e.,  we  may  assert  with  95%  confidence  that  the  mean  value 
lies  between  2.3600  and  3.4000  st./grid  square. 


Save  results  to  a file?  Y/N:  y 
Output  to  file;  BootOut.txt?  Y/N:  n 
Enter  output  file  name:  example!. out 
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Fig.  X2.3  Output  file  for  the  bootstrapping  program  (PC  version) 


EXAMPLES. OUT  - 6/13/95 

Input  data  set  = EXAMPLE3.DAT 
Number  of  grid  openings  = 50 

Example  data  set 


st./grid  square 
For  input  data  set: 


Mean  value 

=s 

2.8800 

bootstrap  data  sets: 

Grand  mean  value 

= 

2.8765 

Lower  95%  confidence 

limit  = 

2.3600 

Upper  95%  confidence 

limit  = 

3.4000 

The  following  would  be  an  appropriate  way  to  state  the  results  of 
this  calculation: 

The  mean  value  for  this  data  set  is  2.8800  st./grid  square.  The  interval 
from  2.3600  to  3.4000  St./grid  square  is  a 95%  confidence  interval  for  t 
mean  value;  i.e.,  we  may  assert  with  95%  confidence  that  the  mean  value 
lies  between  2.3600  and  3.4000  st./grid  square. 


X3.  GENERAL  OUTLINE  OF  THE  BOOTSTRAPPING  PROGRAM 


1.  Read  the  input  data  set 

2.  Determine  the  niunber  of  values  (n)  in  the  input  data  set 

3.  Calculate  the  mean  value  for  the  data  in  the  input  data  set 

4.  LxK)p  for  6000  iterations 

- Loop  for  n times 

- Randomly  choose  a value  from  the  input  data  set 
— Put  the  value  in  a temporary  bootstrap  data  set 

- Replace  the  value  in  the  input  data  set 
— End  loop 

- Determine  the  mean  value  of  the  temporary  bootstrap  data  set 

- Put  the  mean  value  in  a bootstrap  array 

- End  loop 

5.  Sort  the  mean  values  in  the  bootstrap  array 

6.  Choose  values  that  are  150  from  the  top  and  bottom  of  the  bootstrap  array  for  the  lower  and  upper  95% 

confidence  limits 

7.  Determine  the  grand  mean  of  the  values  in  the  bootstrap  array 

8.  Report  the  mean  value  for  the  input  data  set,  the  grand  mean  of  the  bootstrap  values  and  the  lower  and 

upper  95%  confidence  limits. 
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X4.  PROGRAM  FOR  BOOTSTRAPPING  ON  AN  APPLE  COMPUTER 


PROGRAM  BOOT.STRAP; 

USES  Types,  QuickDraw.  Packages,  CursorCtl.  ToolUtils, 

SANE.  PasLibIntf.  StandardFile,  Packinif,  OSUtils; 

CONST 

np  = 6000: 

TYPE 

RealAriayNP  = ARRAY  [l..npj  OF  real; 

IntArravSO  = array  [l-.SO]  of  reai; 
fnamestr  = STRING[631; 


VAR 

out_data 
Ran2Iy 
Ran2Ir 
in.daca 
temp_daia 
n_dat,  ii,  jj  : 
i_ran 
r_ran 

head_x 
strl 

min_val. 
mean_val  : 
mean_2000 


: ReaiArrayNP; 

; longint: 

: ARRAY  [1..971  OF  longint; 

: IntAiraySO:  {structure  values} 

: IntArrayiO; 
integer: 

: longint: 

: real: 

: ARRAY  [1..3]  OF  str255;  (3 

: ARRAY  [1..50)  OF  STRING[10]; 
max_val  : IntArxaySO: 

IntArraySO: 

: IntArraySO: 


median_val 
ji,  n_p 
fhdrinfo 
errCode 
runs,  ir 
tt 

ave_tt 

units 

f_name 

datelnfo 

isec 


: IntArraySO: 
; integer; 
FInfo; 

; integer: 
integer; 

reai: 

: real: 
string[20]; 

: fnamestr: 
DateTimeRec; 
longint: 


lines  of  header  information} 


LABEL 

100: 

FUNCTION  ItsNaN(theString:  str255):  boolean: 

{function  is  a test  to  determine  if  ASCII  characters  are  numbers} 


BEGIN 

ItsNaN  :=  false; 

IF  (ord(ClassExtended(Str2Num(theString))))  < 3 THEN  ItsNaN  = true- 
END: 


FUNCTION  read_file:  integer: 
{function  reads  the  input  data  file} 


VAR 

myReply 

rayTypes 

myErr 


: SFRepiy: 

: SFTypeList: 
: integer: 
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str  : Str255; 

theFile  : text: 

temp  ; extended; 

ij,nt  ; integer; 

f_name  ; fnamestr; 

p : Point; 

cfa  ;CHAR; 

strx  : Stringil]; 

er_NaN  ; boolean: 


BEGIN 


myTypes(0]  :=  TEXT'; 
ch  :=  Chr(ord(9)); 
strx  ;=  ch; 

(♦myForm.Styie  :=  FixedDecimal; 
rnyForm-Digits  ;=  4;») 
p.h  :=  100;  p.v  :=  80; 

SFGetFile(p,  'Select  File  '.  NIL.  1,  myTypes,  NIL.  ray  Reply); 

if  myRepiy.good  then 

BEGIN 

myerr  ;=  0; 

f_name  ;=  myReply.fname; 
reset(theFile,f_name); 
readln(theFile,head_x[l  ]); 
readIn(theFile.head_x[2]) 
readIn(theFile,head_x(3]) 
readln(theFile. units); 

(•readln(theFile.str); 
temp  ;=  Str2Num(str); 
n_ciat  :=  Num2Integer(temp); 

FOR  ij  :=  1 to  n_dat  do  BEGIN*) 
ij;=0; 
tt  :=  0; 

REPEAT 


[read  three  header  lines} 


{read  the  units  for  the  values} 


readln(theFile,str); 

(read  in  the  identification  for  the  grid  openings  analyzed 
[maximum  eight  charactgers]  and  the  value  for  the  number  or 
concentration  of  asbestos  structures} 
ij  :=ij+  1; 

nt  :=  Pos(strx.  str)  -1; 

sixUij]  ;=  Copy(str,  1,  nt); 

nt  :=  nt  + I ; 

Delete  (str,  1,  nt); 

er_NaN  :=  ItsNaN(str); 

if  er_NaN  = false  then  begin 
temp  ;=  Str2Num(str); 
in_data[ij]  :=  temp; 

tt  :=  tt  + in_data[ij];  (sum  of  data  values} 

(*Num2Dec(myForm,  temp,  in_data[ij]);*) 

end 

else  myErr  ;=  2; 

UNTIL  (EOF(theFile)  = true)  OR  (myerr  = 2); 

n_dat  :=  ij;  (number  of  values  in  the  input  data  set} 

ave_tt  :=  tt/ij;  (calculate  the  mean  value  for  data  in  the 
input  data  set} 
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ciose(theFile); 

{myErr  :=  0;) 

END 

ELSE 

myErr  :=  1; 
read_file  :=  myErr; 

END. 

PROCEDURE  write_xfile;  {procedure  for  writing  the  output  file} 

VAR 

my_Reply  ; SFReply; 
myErr  : integer; 

f_nam  : fnamestr; 

the_File  : text; 

temp  : extended; 

Ut  vol_num  : integer; 

my  Form  : DecForm; 

str,  str2  : DecStr, 

str3.  str4  : DecStn 

str5  : DecStr. 

str6  : Str255; 

p Point; 

BEGIN 

my_Reply.fType  :=  TEXT'; 
p.fa  :=  100;  p.v  ;=  80; 
getdatetime(isec); 

secs2Date(isec.  datelnfo); 

f_nam  :=  ’Bootruns.txt'; 

SFPutFile(  p/Save  an  Output  File',  f_nam,  NIL,  my_Reply); 
if  my_Reply.Good  then  BEGIN 
f_nam  :=  my_Reply.fname; 
vol_num  :=  my_Reply.VRefNum; 
rewrite(the_File,f_nam); 

{ fndrlnfo.fdcreator  ;=  'XCEL';  } 

fndrlnfo.fdcreator  :=  'MSWD'; 
fndrinfo.fdtype  :=  TEXT; 
fndrlnfo.fdlhags  :=  0{$100}; 
myErr  :=  setfinfo(f_nam,  vol_num,  fndrinfo); 

writeln(the_File,'  ',f_nani,’  - ’,dateInfo.month;2,'/',dateInfo.day;2,'/’,dateInfo.year- 

1900;2  ); 

writeln(the_File); 

writein(the_File, 'Input  data  set  ='  ',chr(9),f_name); 

writeln(the_File,’Number  of  grid  openings  = ',n_dat); 

writeIn(the_File); 

writeln(the_File,head_x[  I ]); 

writeln(the_File,head_x[2]); 

writeIn(the_File,head_x[3]); 

writeln(the_File, units); 

str6  :=  concatCMean  value',  chr(9),  'Grand  Mean',  chr(9),  'Median  value’,  chr(9),  'Lower 
limit'.  chr(9).  'Upper  limit'); 

{ writeln(the_Fi]e.str6);  } 

myForm.Style  :=  FixedDecimal; 
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myForm.Digits  :=  3; 

FOR  ij  :=  1 to  runs  DO  begin 

Num2str(myForm.  mean_val[ij],  str); 

Num2str(myForin,  mean_2000[ij],  str2); 

Num2str(myFonn.  median_val[ij],  str3); 

Num2str(myFonn.  min_val[ij],  str4); 

Num2str(myFonn.  max_val[ij],  str5); 

Sff6  :=  concat(str.  chr(9),  str2,  chr(9),  str3,  chr(9),  str4,  chr(9),  str5); 
{ writelnCthe.File.strb);  } 

writeln(the_File); 

writeln  (the_File/  For  input  data  set:  ' 

writeln  (the_FiIe,'  Mean  value  = ',chr(9),meaii_vairij]:9:4)- 

writeIn(the_FiIe); 

wntein  (the_File.’  For  bootstrap  data  sets:  '); 

writeln  (the_File.’  Grand  mean  value  = ',chr(9),mean_2000[ij]:9:4); 

wntein  (the_File.  Lower  95%  confidence  limit  = ',min_val[ij]:9:4); 

wntein  (the_File.'  Upper  95%  confidence  limit  = ',inax_val[ij]:9:4); 

writeln(the_File); 

writeln  (the_File.  The  following  would  be  an  appropriate  way  to  state  the 
'results  of); 

wntein  (the_File.'  this  calculation:  ’); 
writein(the_File); 

wntein(the_File,’  The  mean  value  for  this  data  set  is 
mean_val[ij]:9:4.’  units.'.  The  interval'); 
writeln  (the_File.'  from'.min_val(ij]:9:4, 

to  ',max_val[ij]:9:4.’  '.units.'  is  a 95%  confidence' 

,'  interval  for  the'); 

writeln  (the_File,'  mean  value;  i.e..  we  may' 

, assert  with  95%  confidence  that  the  true  mean  value  '); 
writeln  (the_FiIe,'  lies  'oetween  ',min_val[ij]:9:4.’  and  ’,max_val[ij];9;4,’ 
units.'.'); 

writein(the_File); 

B®;  (FOR  ij  ;=  1 to  runs} 
close(the_File); 

niyErr  :=  Flush Vol(NIL.  vol_num); 

END. 

END, 

FUNCTION  ran2(VAR  idum:  longint):  real; 

(function  to  generate  random  numbers  for  the  bootstrapping  program) 

CONST 

m = 714025; 
ia  = 1366; 
ic  = 150889; 
rm  = 1.40051 12e-6; 

VAR 

j:  integer; 

BEGIN 

IF  idum  < 0 THEN  BEGIN 
idum  :=  (ic-idum)  MOD  m; 

FOR  j :=  I TO  97  DO  BEGIN 
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idum  :=  (ia*idum+ic)  MOD  m; 
Ran2Ir(j]  :=  idum: 

END; 

idum  :=  (ia*idum+ic)  MOD  m; 
Ran2Iy  :=  idum; 

END; 

j :=  1 + (97*Ran2Iy)  DIV  m; 

( IF  0 > 97)  OR  (j  < 1)  then  begin 

writelnCpause  in  routine  RAN2'); 
readln; 

END;} 

Ran2Iy  :=  Ran2Irlj]; 
ran2  :=  Ran2Iy*rm: 
idum  :=  (ia*idum+ic)  MOD  m; 
Ran2Ir(j]  ;=  idum; 

END. 


FUNCTION  calc_mean(na  : integer;  z : RealArrayNP):  REAL; 

(function  to  calculate  the  grand  mean  of  the  means  for  ail  bootstrap  data  sets) 

VAR 

s ; real; 

i : integer; 

BEGIN 

s :=  0.0; 

FOR  i :=  1 to  na  DO 
s ;=  s + z(i]; 
calc_mean  ;=  s/na; 

END, 

FUNCTION  caic_mean50(na  : integer;  z ; IntArraySO):  REAL: 

(function  to  calculate  the  mean  value  for  a single  bootstrap  data  set} 

VAR 

s : real; 

i : integer; 

BEGIN 

s ;=  0.0; 

FOR  i :=  1 to  na  DO 
s :=  s + zfi]; 
calc_mean50  :=  s/na; 

esiD. 

PROCEDURE  son_data(n;  integer; 

VAR  ra:  RealArrayNP); 

(procedure  to  sort  the  mean  values  determined  for  the  6000  bootstrap  data  sets} 

LABEL  99; 

VAR 

l,j,ir,i;  integer, 
rra:  real; 

BEGIN 

1 ;=  (n  DIV  2)+I; 
ir  ;=  n; 

WHILE  true  DO  BEGIN 
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IF  1 > I THEN  BEGLN 
1 :=  1-1; 
rra  :=  ra[l] 

END 

ELSE  BEGIN 
rra  :=  ra(ir]; 
ra[ir]  :=  ra(l]; 
ir  ;=  ir-1; 

' IF  ir  = 1 THEN  BEGIN 
ra(l]  :=  rra: 

GOTO  99; 

END; 

END; 

1 :=I; 
j:=  1+1; 

WHILE  j <=  ir  DO  BEGIN 
IFj  <ir  THEN 
IF  ra(j]  < ra(j+I ) THEN 
J :=j+i; 

IF  rra  < ra[j]  THEN  BEGIN 
ra[i]  ;=  ra(j]; 
i :=j; 
j :=j+j; 

END 
ELSE 
j :=  ir+I; 

END; 

ra(i]  :=  rra; 

END; 

99: 

END, 

BEGIN  {Main  program} 

writelnC  NIST  data  analysis  program.'); 

( writelnC  Enter  number  of  iterations  (max  6000,  min  40).’); 
readln(n_p); 

writeinC  Enter  number  of  runs  (max  50).’); 
readln(runs);  } 

InitCursorCtl(NIL); 

n_p  :=  6(XK); 
runs  :=  1; 

writelnC’  Hit  Return  and  select  a file  to  open.'); 
readln; 

100: 

errCode  :=  read_file; 
if  (errCode  = 0)  then  begin 

writelnC’  Calculating  bootstrap  - ’); 
SetCursorCGetCursorCWATCHCURSOR)''''); 
i_ran  :=  -1  (TickCountj; 

(This  line  uses  the  tick  count  to  generate  an  initial  random  starting  value 
for  the  random  number  generator} 

if  Cn_p  < 40)  then  n_p  :=  40; 

FOR  ir  ;=  1 to  runs  DO  BEGIN 
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FOR  ii  ;=  I TO  n_p  DO 
BEGIN 

RoiaieCursor(TickCount); 

FOR  ji  ;=  1 TO  n_dat  DO 
BEGIN 

r_ran  :=  ran2(i_ran); 
jj  ;=  TRUNC(r_ran  * n_dat)  + 1 ; 
if  jj  > n_dat  then  jj  :=  n_dat; 
teinp_data[ji]  :=  in_data|jj]; 

END, 

out_data(ii]  :=  caic_mean50(n_dat.temp_data); 
(FOR  ii  :=  1 TO  n_p} 

mean_vai[ir]  ;=  calc_mean50(n_dat.in_data); 

mean_2000[ir]  :=  calc_mean(n_p,out_data); 

sort_data(n_p,oui_data); 

ji  :=  TRUNC(n_p  * 0.025); 

min_val[irl  :=  out_data[ji]; 

max_val[ir]  ;=  out_data[n_p  - ji]; 

ii  :=  n_p  DIV  2; 

median_val[ir]  :=  out_data[ii]; 

InitCursor; 

writein; 

{This  information  is  given  on  the  screen  as  the  program  is  run} 
writein: 

writein  ('  For  input  data  set:  ’); 

writein  ('  Mean  value  = ',mean_val(ir]:9:4); 

writein; 

writein  (’  For  bootstrap  data  sets:  '); 

writein  (’  Grand  mean  value  = ■,mean_2000[ir]:9:4); 

writein  (’  Lower  95%  confidence  limit  = ’,niin_val[irl:9:4); 

writein  (’  Upper  95%  confidence  limit  = ’,max_val[ir]:9:4); 

writein: 
writein; 

writein  ('  The  following  would  be  an  appropriate  way  to  state  the 
'results  of); 

writein  ('  this  calculation;  '); 
writein; 

writelnC  The  mean  value  for  this  data  set  is 

mean_val(ir]:9:4,’  units,'.  The  interval'); 
writein  ('  from’,min_val[ir]:9;4, 

' to  ',max_val[ir]:9:4,’  '.units,'  is  a 95%  confidence' 

,'  interval  for  the'); 
writein  ('  mean  value;  i.e.,  we  may' 

,'  assen  with  95%  confidence  that  the  true  mean  value  '); 
writein  ('  lies  between',min_val[ir]:9:4.'  and’,max_val[ir]:9:4,'  ', 
units,'.'); 

END, 

writein; 

writelnC  Hit  Return  and  select  a file  for  output  or  cancel.'); 
readln; 

write_xfile;  (write  the  output  file} 

end 

else  if  (errCode  = 2)  then  goto  100; 
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X5.  PROGRAM  FOR  BOOTSTRAPPING  ON  A PC 


uses  Dos,  Crt,  windos,  strings,  printer; 
const 


np  = 6000; 


TYPE  RealArrayNP  = ARRAY  [1  ..np]  OF  real; 
IntAiraySO  = ARRAY  [1  ..50]  of  real; 
fiiameStr  = ARRAY  [0..62]  of  char; 


VAR 


opti  : boolean; 

optr  : boolean; 


def_prt  : boolean; 

out_data  : RealArrayNP; 

Ran2Iy  : longint; 

Ranllr  : ARRAY  [1  ..97]  OF  longint; 

in_data  : IntArraySO;  {structure  values} 

temp_data  : IntArraySO; 

n_dat,  ii,  jj  : integer, 

i_ran  : longint; 

r_ran  : real; 

head_x  ; ARRAY  [1..3]  OF  STRING[255];  {3  lines  of  header  information} 

strl  :ARRAY[1..50]OFSTRING[10]; 

min_val,  max_val  : IntArraySO; 

mean_val  ; IntArraySO; 

mean_2000  ; IntArraySO; 

median_val  : IntArraySO; 

ji,  n_p  : integer, 

errCode  : integer, 

code  : integer; 

runs,  ir  : integer; 

f_name  : fiiameStr; 

stra  ; string[I5]; 

strb  :string[15]; 

strc  :string[15]; 

strd  :string[15]; 

stre  ; string[15]; 

strf  ; string[60]; 

ave_tt  : real; 
the_File  : text; 

f_nam  ; fiiameStr; 

1st  ; text ; 

prt  ; Boolean; 

units  : string[20]; 


FUNCTION  read_file;  integer; 

{this  function  reads  the  input  data  file} 


VAR 


myErr 

stm 


: integer, 

: STRING[8]; 
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theFile 

: text; 

temp 

: real; 

ij,nt 

: integer. 

ch 

:CHAR; 

tt 

: real; 

tl 

:string[12];  {sum  of  data} 

11 

: integer; 

S;  array [0..fsPathName]  of  Char;  {used  for  fileSearch} 

label  10; 
label  20; 

BEGIN 

ch  :=  Chr(ord(9)); 

10:  FileSearch(S,  f_name,");{  check  to  see  if  file  exists} 

ifS[0]  = 

#0  then 

Begin 

{file  does  not  exist} 

Writeln(f_name,'  not  found  Please  input  file  name: '); 
readln(f_name); 
goto  10; 
end; 

assign  (theFile/_name); 
reset(theFile);  {open  the  file  for  reading} 

20:  readln(theFile^ead_x[l  ]);  {read  the  three  header  lines} 
readln(theFile4iead_x[2]); 
readln(theF  ile  Jiead_x  [3  ] ) ; 
readln(theFile,units); 

tt:=0; 
ij  :=  0; 

REPEAT 
ij  :=ij  + l; 

tt  :=tt  + m_data[ij];  {sum  of  data} 
readln(theFile,  stm,in_data[ij  ] ) ; 

{read  in  the  identification  for  the  grid  openings  analyzed 
[maximum  eight  characters]  and  the  value  for  the  number 
or  concentration  of  asbestos  structures} 

tt  :=  tt  + in_data[ij];  {sum  of  data  values} 

UNTIL  EOF(theFile)  = true; 

ave_tt  :=tt/ij;  {calculate  the  mean  value  for  data 
in  the  input  data  set} 

n_dat  :=  ij;  {number  of  values  in  the  input  data  set} 

close(theFile); 
myErr  :=  0; 
read_file  :=  myErr, 

END; 
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PROCEDURE  •write_xfile;  {procedure  for  writing  the  output  file} 
VAR 

myErr  : integer; 

f_nam  : fiiameStr, 

temp  : real; 
ij,  vol_num  : integer; 
str3  : String[255]; 

S : ARRAY[0..fsPathName]  of  Char; 

yr,  mon,  day,  dow  : word; 
const 

def_file  : ARRAY[0..  1 2]  of  char  = 'BootOut.txt'; 

label  10; 
label  999; 

BEGIN 

ij;=l; 

if  (def_prt)  then 
begin 

{The  information  below  is  written  to  the  screen  as  the  program  is  run} 
writeln; 

writeln  ('  For  input  data  set: '); 

writeln  ('  Mean  value  = ',mean_val[ij]:9:4); 

writeln; 

writeln  ('  For  bootstrap  data  sets: '); 

writeln  ('  Grand  mean  value  = '4nean_2000[ij]:9:4); 

writeln  ('  Lower  95%  confidence  limit  = ',min_val[ij]:9:4); 

writeln  ('  Upper  95%  confidence  limit  = ',max_val[ij]:9:4); 

writeln; 

writeln; 

writeln  ('  The  following  would  be  an  appropriate  way  to  state  the 
'results  of); 

writeln  ('  this  calculation: '); 
writeln; 

writelnC  The  mean  value  for  this  data  set  is 
mean_val[ij]:9:4,' units,'.  The  interval'); 
writeln  ('  fi‘om',min_val[ij]:9:4, 

' to  ',max_val[ij]:9:4,'  ',units,'  is  a 95%  confidence' 

,'  interval  for  the'); 
writeln  ('  mean  value;  i.e.,  we  may' 

,'  assert  with  95%  confidence  that  the  true  mean  value '); 
writeln  ('  hes  between',min_val[ij]:9:4,'  and',max_val[ij]:9:4,' ', 
units,'.'); 
end; 
writeln; 
writeln; 

write('  Save  results  to  a file?  Y/N: '); 

readln(str3[l]); 
if  (pit)  then 
begin 

writeln(lst); 

writeln(lst); 
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wnteln  (1st,'  The  following  would  be  an  appropriate  way  to  state  the 
'results  of); 

writeln  (1st,'  this  calculation: '); 
writeln(lst); 

writeln(lst,'  The  mean  value  for  this  data  set  is ', 
mean_val[ij]:9;4,' ',  units,'.  The  interval'); 
writeln  Ost,'  from',min_val[ij]:9;4, 

' to  ',max_val[ij]:9:4,'  ',units,'  is  a 95%  confidence' 

,'  interval  for  the'); 

writeln  (1st,'  mean  value;  i.e.,  we  may' 

,'  assert  with  95%  confidence  that  the  true  mean  value '); 
writeln  (1st,'  lies  between  'jnin_val[ij]:9:4,'  and  ',max_val[ij];9:4,' ', 
units,'.'); 
writeln(lst); 
writeln(lst); 

writeln(lst,'  Save  results  to  a file?  Y/N:  ',str3[l]); 
end; 

if(UpCase(str3[l])  = "N')  then  goto  999; 
writeC  Output  to  file:  BootOut.M?  Y/N '); 
readln(str3[l]); 

if  (pit)  then  writeln(lst,'  Output  to  file:  BootOut.txt?  Y/N:  ',str3[l]); 
10:  if(UpCase(str3[l])  = 'Y')  then  {use  BootOut.txt} 
begin  {check  to  see  if  file  already  exists} 

FileSearch(S,def_file,  GetEnvVar(TATH')); 
if  S[0]  o #0  then 
Begin 

writeC  ',def_file,'  already  exists.  Do  you  want  to  replace  it?  Y/N: '); 
readln(str3[l]); 

if  (pit)  then  writeln(lst,'  ',def_file,'  already  exists.', 

' Do  you  want  to  replace  it?  Y/N:  ',str3[l]); 
if  (UpCase(str3[l])  = 'N')  then  goto  10; 
writelnC  Output  sent  to  ',def_file); 
if(prt)  then  writelnGst,'  Output  sent  to  ',def_file); 
assign(the_F  ile,def_file) ; 
end; 


end 

else 

begin 

writeC  Enter  output  file  name: '); 
readln(f_nam); 

FileSearch(ST_nam,  GetEnvVar(PATH')); 

if  (prt)  then  write(lst,'  Enter  output  file  name:  '/_nam); 
if  S[0]  o #0  then 
Begin 

write("T_nam,' already  exists.  Do  you  want  to  replace  it?  Y/N'); 
read[ln(str3[l]); 

if  (prt)  then  write('  'jf_nam,'  already  exists.  Do  you  want  to  replace  it?  Y/N:  ',str3[l]) 
if  (UpCase(str3[l  ])  = 'N')  then  goto  10; 
assign(the_F  ile,def_file) ; 
end; 

assign(the_F  ileT_nam) ; 
end; 
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getdate(yr,  mon,  day,  dow);  {get  date  for  output  file} 

rewrite(the_File); 

wTiteln(the_File); 

if  (UpCase(str3  [ 1 ])  = 'Y')  then  writeln(the_File,' 

StrUpper(def_file),'  - ',mon,'/',  day,'/*,  yr- 1900) 

else 

writeln(the_File,' 

StrUpper(f_nam), ' - ',mon,'/’,day,'/',yr-1900); 
writeln(the_File); 

writeln(the_File,'  Input  data  set  = ',StrUpper(f_name)); 
writeln(the_File,'  Number  of  grid  openings  = 'm_dat); 

if(  runs  > 1 ) then 
begin 

writeln(the_File,'  Number  of  iterations  (bootstrap  data  sets)  = ',njp); 

writeln(the_File,'  Number  of  runs  = ' juns); 

end; 

writeln(the_File); 
writeln(the_File,'  'diead_x[l]); 
writeln(the_File,’  '4iead_x[2]); 
vmteln(the_File,'  ’diead_x[3]); 
writeln(the_File,'  ',umts); 
vmteln(the_File); 

if((runs=  1)  and  (n_p  = 6000))  then 
begin 

writeln  (the_File,'  For  input  data  set: '); 

writeln  (the_File,'  Mean  value  = 'miean_val[ij]:9:4); 

writeln(the_File); 

writeln  (the_File,'  For  bootstrap  data  sets: '); 
writeln  (the_File,'  Grand  mean  value  = ',mean_2000[ij]:9:4); 
writeln  (the_File,'  Lower  95%  confidence  limit  = ',min_val[ij]:9:4); 
writeln  (the_File,'  Upper  95%  confidence  limit  = ',max_val[ij]:9:4); 

writeln(the_File); 

writeln(the_File); 

writeln  (the_File,'  The  following  would  be  an  appropriate  way  to  state  the 
'results  of); 

writeln  (the_File,'  this  calculation: '); 
writeln(the_File); 

writeln(the_File,'  The  mean  value  for  this  data  set  is 
mean_val[ij]:9:4,' units,'.  The  interval'); 
writeln  (the_File,'  fi-om  ',min_val[ij]:9:4, 

' to  ',max_val[ij]:9:4,'  ',units,'  is  a 95%  confidence' 

,'  interval  for  the'); 

writeln  (the_File,'  mean  value;  i.e.,  we  may' 

,'  assert  with  95%  bonfidence  that  the  true  mean  value '); 
writeln  (the_File,'  lies  between  'miin_val[ij]:9:4,'  and  ',max_val[ij]:9:4,' ', 
units,'.'); 

writeln(the_File); 
writeln(the_F  ile); 
end; 

str3  :=  concatCMean  value',  chr(9),  'Grand  Mean',  chr(9). 

Median  value',  chr(9).  Tower  limit',chr(9),  'Upper  limit'); 
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writeln; 
if  ( opti)  then 
begin 

writeln(the_File); 

writeln(the_File,'Run  Mean  Value  Grand  Mean 
."Median  Lower  Upper'); 

writeln(the_FiIe,'  (of  input  bootstrap  Limit 

writeln(the_File,'  counts)  data  sets'); 

writeln(the_File); 

FOR  ij  :=  1 to  runs  DO  begin 

str(mean_val[ij] : 1 0:4,stra); 
str(mean_2000[ij] : 1 0:4,strb)  ; 
str(median_val[ij] ; 1 0:4,strc); 
str(min_val[ij] : 1 0:4,strd); 
str(max_val[ij] : 1 0:4,stre); 

strf  :=  concat(stra,chr(9),strb,chr(9),strc,chr(9),strd,chr(9),stre); 
writeln  (the_file,'  ',ij,'  '.strf); 

end; 


end 

else 

begin 

end; 

close(the_File); 

999:  END; 


FUNCTION  ran2(VAR  idum:  longint):  real; 

{function  to  generate  random  numbers  for  the  bootstrapping  program} 

CONST  m =714025; 
ia=  1366; 
ic=  150889; 
rm=  1.40051 12e-6; 

VAR 

j:  integer; 

BEGIN 

IF  idum  < 0 THEN 
BEGIN 

idum  :=  (ic-idum)  MOD  m; 

FOR  j :=  1 TO  97 
DO  BEGIN 

idum  :=  (ia*idum+ic)  MOD  m; 

Ran2Ir[j]  :=  idum; 

END; 

idum  :=  (ia*idum+ic)  MOD  m; 

Ran2Iy  :=  idum; 

END; 
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Limit') 


j :=  1 + (97*Ran2Iy)  DIV  m; 

Ran2Iy  :=  Ran2Ir[j]; 

ran2  ;=  Ran2Iy*nn; 

idum  :=  (ia*idum+ic)  MOD  m; 

Ran2Ir[j]  :=  idum; 

END; 

FUNCTION  calc_mean(na  : integer;  z : RealArrayNP):  REAL; 
{function  to  calculate  the  grand  mean  for  all  the  bootstrap  data  sets} 

VAR 

s : real; 

i : integer, 

BEGIN 

s:=0.0; 

FOR  i :=  1 to  na  DO 
s :=  s + z[i]; 
calc_mean  :=  s/na; 

END; 

FUNCTION  calc_mean50(na  : integer,  z : IntArraySO):  REAL; 
{flmction  to  calculate  the  mean  value  for  a single  bootstrp  data  set) 

VAR 

s : real; 

i : integer; 

BEGIN 
s :=  0.0; 

FOR  i ;=  1 to  na  DO 
s :=  s + z[i]; 
calc_mean50  :=  s/na; 

END; 


PROCEDURE  sort_data(n:  integer; 

VAR  ra:  RealArrayNP); 

{procedure  to  sort  the  mean  values  determined  for  the  6000  bootstrap 
data  sets) 

LABEL  99; 

VAR 

lj,ir,i:  integer; 
ira:  real; 

BEGIN 

I:=(nDIV  2)+l; 
ir  :=  n; 

WHILE  true  DO  BEGIN 
IF1>1  THEN  BEGIN 
1:=1-1; 
rra  :=  ra[l] 

END 

ELSE  BEGIN 
rra  :=ra[ir]; 
ra[ir]  :=ra[l]; 
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ir  :=  ir-l; 

IFir=  1 THENBEGIN 
ra[l]  :=  rra; 

GOTO  99; 

END; 

END; 
i:=  1; 
j:=l+l; 

WHILE  j‘  <=  ir  DO  BEGIN 
IF  j < ir  THEN 
IF  ra[j]  < ra[j+l]  THEN 

j :=j+i; 

IF  rra  < ra[j]  THEN  BEGIN 
ra[i]  :=  ra[j]; 

i:=j; 
j ~j+j; 

END 
ELSE 
j ;=ir+l; 

END; 
ra[i]  ;=  rra; 

END; 

99; 

END; 


var 

hr,  min,  sec,  hund  : word; 
opts  : string; 

label  10; 

BEGIN  {Main  program} 

{note:  opti,  optr,  prt  refer  to  additional  features  not  described  in 
thisNISTIR} 

opti  :=  False; 
optr  :=  False; 
prt  :=  False; 
def_prt  :=  True; 
n_p  ;=6000; 
runs  :=  1; 

clrScr, 

for  ii  :=  1 to  ParamCount 
do  begin 

opts  :=  concat(opts,  ParamStr(ii))  ; 
end; 

ifCCPosCR^opts)  > 0)  or  (Pos('r',opts)  > 0))  then  optr  :=  True; 
if((Pos('r,opts)  > 0)  or  (Pos('i',opts)  > 0))  then  opti  ;=  True; 
if((Pos(P',opts)  > 0)  or  (Pos('p',opts)  > 0))  then  prt  ;=  True; 
if  (optr  or  opti)  then  def_prt  :=  False; 
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if  (prt)  then 
begin 

writelnC  Print  the  screen  inhouse  option'); 

writelnC  Make  sure  Printer  is  online  Hit  return  when  ready'); 

readln; 

assign(Lst,'lptr); 

rewrite(lst); 

end; 

writelnC  NIST  data  analysis  program'); 
writeln; 

writeC  Input  the  file  name: '); 
readln(f_name); 

if  (opti)  then 
begin 

writeC  Enter  number  of  iterations  (max  6000).: '); 

readln(n_p); 

end; 

if  (optr)  then 
begin 

writeC  Enter  number  of  runs  (max  50).: '); 
readln(runs); 
end; 

if  prt  then 
begin 
{$!-} 

10:  writeln(lst,'  NIST  data  analysis  program'); 

{$!+} 

if  lOResult  = 1 60  then 
begin 

writelnC  Printer  is  off-line?'); 
writeln  ('  hit  Return  when  ready'); 
readln; 
goto  10; 
end; 

writeln(lst); 

write(lst,'  Input  the  file  name: '); 

writeln(lstT_name); 

writelnOst); 

if  (opti)  then 
begin 

write(lst,'  Enter  number  of  iterations  (max  6000).: '); 
writeln(lsUi_p); 

write(lst,'  Enter  number  of  runs  (max  50).: '); 
writelnOstjuns); 
end; 
end; 

clrScr; 

writelnC  Calculating  bootstrap  - '); 
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if  (pit)  then  writeln(lst,'  Calculating  bootstrap  - '); 
errCode  :=  read_file; 

if  (errCode  = 0)  then  begin  {file  OK} 

{The  next  three  lines  use  the  time  to  generate  an  initial 
random  starting  value  for  the  random  number  generator} 

i_ran  :=  -1; 

GetT  ime(hr4nin,secdiund); 
i_ran  :=  hund  * min  * sec; 

FOR  ir  :=  1 to  runs  DO  BEGIN 
FOR  ii  :=  1 TO  n_p  DO  BEGIN 

FOR  ji  :=  1 TO  n_dat  DO  BEGIN 
r_ran  :=  ran2(i_ran); 
jj  :=  TRlJNC(r_ran  * n_dat)  + 1 ; 
if  jj  > n_dat  then  jj  :=  n_dat; 
temp_data[ji]  :=  in_dataljj]; 

END; 

out_data[ii]  :=  calc_mean50(n_dat,temp_data); 

END; 

mean_val[ir]  :=  calc_mean50(n_dat,in_data); 

mean_2000[ir]  :=  calc_mean(n_p,  out_data); 

sort_data(n_p,out_data); 

ji  :=  TRUNC(n_p  * 0.025); 

min_val[ir]  :=  out_data[ji]; 

max_val[ir]  :=  out_data[nj)  - ji]; 

ii:=n_pDIV  2; 

median_val[ir]  :=out_data[ii]; 
str(mean_val[ir]:  10:4,stra); 
str(mean_2000[ir]:10:4,strb)  ; 
str(min_val[ir] : 1 0:4,strd); 
str(max_val[ir] ; 1 0;4,stre); 

strf  :=  concat(stra,chr(9),strb,chr(9),strc,chr(9),strd,chr(9),stre); 

if  (optr)  then 
Begin 
writeln; 

writelnC  Run  number:  ',ir); 
writelnC  Mean  value  is:  ',chr(9),stra); 
writelnC  Mean  value  of  ',n_p,'  is:  ',strb); 
writelnC  Lower  limit  is:  ',chr(9),strd); 

writelnC  Upper  limit  is:  ',chr(9),stre); 

end; 

if  (pit)  then 
Begin 

if  (opti)  then 

Begin 

writeln(lst); 

writeln(lst,'  Run  number:  ',ir); 
writeln(lst,'  Mean  value  is:  ',chr(9),stra); 
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writelnOst,'  Mean  value  of is:  ',chr(9),strb); 
writeln(lst,'  Median  value  is:  ',chr(9),strc); 

vvriteln(lst,'  Lower  limit  is:  ',chr(9),strd); 

writeln(lst,'  Upper  limit  is:  ',chr(9),stre); 

end; 

if  (prt  and  def_prt)  then 
begin 

writeln(lst); 

writeln  (1st,'  For  input  data  set: '); 
writeln  (1st,'  Mean  value  = ', 

mean_val[ir]  :9 :4); 
writeln(lst); 

writeln  (1st,'  For  bootstrap  data  sets: '); 
writeln  (1st,'  Grand  mean  value  = ', 
mean_2000[ir]  :9:4); 

writeln  (1st,'  Lower  95%  confidence  limit  = ', 
mm_val[ir]:9:4); 

writeln  (1st,'  Upper  95%  confidence  limit  = ', 
max_val[ir]:9:4); 

end; 

end; 

END; 

write_xfile;  { write  the  output  file } 
end; 

if  prt  then  writeln(lst,chr(12)); 
if  prt  then  close(lst); 

END. 
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